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iH : Abstract 



The propagation of bound optical waves along the surface of a one-dimensional (1-D) photonic 
crystal (PC) structure is considered. A unified description of the waves in 1-D PCs for both s- 



^ ■ and p-polarizations is done via an impedance approach. A general dispersion relation that is valid 



for optical surface waves with both polarizations is obtained, and conditions are presented for 
long-range propagation of plasmon-polariton waves in nanofilms (including lossy ones) deposited 
on the top of the 1-D PC structure. A method is described for designing 1-D PC structures to 
fulfill the conditions required for the existence of the surface mode with a particular wavevector 
at a particular wavelength. It is shown that the propagation length of the long-range surface 
plasmon-polaritons can be maximized by wavelength tuning, which introduces a slight asymmetry 



>, 
^^ . in the system 

o 



PACS numbers: 42.70.Qs, 73.20.Mf, 78.67.-n, 78.68. +m 



I. INTRODUCTION 

Optical surface waves (SWs) are excitations of electromagnetic (EM) modes that are 
bound to the interface between two media. The maximum EM field strength of the SWs is 
located near the interface, so these waves have a high sensitivity to the surface condition. 
The SWs that have been employed in the widest array of applications are surface plasmon 
polariton (SPP) waves, which are p-polarized optical surface waves propagated along a 
metal-dielectric interface^. The sensitivity of SPPs has been used in many surface plasmon 
resonance (SPR) applications in which a shift of an SPR dip is measured, from biosensors 
used for detecting biomolecules in a liquid, to gas sensors for detecting trace impurities in 
the air-. 

A limiting factor for SPR sensitivity in these applications is the limited propagation 
length of SPPs due to a strong intrinsic damping of their EM field in metal. Even when 
the "best plasmonic" metals, such as silver and gold, are used, the SPP propagation length 
is only about ten micrometers in the optical frequency range. Other metals do not prac- 
tically support any SPP propagation at visible frequencies. One way to increase the SPP 
propagation length and, consequently, the SPR sensitivity is to use long-range SPPs (LR- 
SPPs), which can be achieved using a thin metal film embedded between two dielectrics 
with identical refractive indices (RIs)^"-. 

In the work-, another method for the excitation of LRSPPs was reported, in which the 
thin metal film was embedded between the medium being studied (with any RI, even a gas 
medium) and the 1-D photonic crystal. Photonic crystals (PCs) are materials that possess a 
periodic modulation of their refraction indices on the scale of the wavelength of light^. Such 
materials can exhibit photonic band gaps that are very much like the electronic band gaps 
for electron waves traveling in the periodic potential of the crystal. In both cases, frequency 
intervals exist in which wave propagation is forbidden. This analogy may be extended^ to 
include surface levels, which can exist in band gaps of electronic crystals. In PCs, they 
correspond to optical surface waves with dispersion curves located inside the photonic band 
gap. 

The one-dimensional photonic crystal (1-D PC) is a simple periodic multi-layer stack. 
Optical surface modes in 1-D PCs were studied in the 1970s, both theoretically^ and experi- 
mentallyJ^. Twenty years later, the excitation of optical surface waves in a Kretschmann-like 



configuration was demonstratedii. A scheme of tlie Kretscliniann-like excitation of PC SWs 
is presented in Fig. [TJ In recent years, tlie PC SWs liave been used in ever-widening applica- 
tions in tlie field of optical sensorsi^i"— . In contrast to SPPs, both p-polarized and s-polarized 
optical surface waves (with a dielectric final layer of the 1-D PC)^^^ can be used in PC SW 
sensor applications. However, in applications in which the LRSPP propagates in a metal 
nanofilm deposited on an appropriate 1-D PC surface (e.g., in hydrogen detection^^ii^), only 
p-polarized waves can be used. Therefore, a unified theoretical description of PC SWs for 
both polarizations would be very useful, and such a description is presented in this article. 
Additionally, the current contribution describes conditions for the propagation of LRSPP in 
metal nanofilms (including lossy ones) that are deposited on the top of the 1-D PC structure 
and, thereby, provides the theoretical background for the works^ii^ii^. 
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FIG. 1: Excitation of optical surface waves in a Kretschmann-like scheme. 



II. DISPERSION RELATION FOR S- AND P-POLARIZED SURFACE WAVES 
IN IMPEDANCE TERMS 

A. Impedance approach: basic definitions and recursion relation for a multilayer 

The "characteristic impedance" of an optical medium is the ratio of the electric field 
amplitude to the magnetic field amplitude in this medium, i.e., ^char = E/H = 1/n. The 
concept of impedance in the optics of homogeneous layers is based on a mathematical analogy 
between a cascade of transmission line sections and a multilayer optical coating. 
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FIG. 2: Reflection and transmission from a single interface. 



For reflection from a plane interface, the useful value is the "normal impedance"— i^^ Z 
which is the ratio of the tangential components of the electric field to the magnetic field: 
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Impedances for the s-polarized wave (in which the electric field vector is orthogonal to the 
incident plane — TE wave) and for the p-polarized wave (in which the electric field vector 
is parallel to the incident plane — TM wave) are correspondingly: 
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For the interface between semi- infinite media and 1 (Figure [2]), Fresnel's formula for 
reflection coefficients has a very simple form in the impedance terms: 

fl = |a^. (4) 

where Zq^ is the normal impedance of medium j, given by (E]) or ([3]). Hereafter, the use of R 
and Z (without subscripts s or p) means that the equation holds for both polarizations when 
the corresponding impedances Z^ or Zp are inserted. Fresnel's formulas for transmission 
coefficients are as follows: 

T, = -2—^f^-- (5) 

Zs{l) + Zs(0) 
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T, = -2^ ^ ^'^'l . (6) 



In Figure [21 the labeling of the s-polarized wave in the first and the second media by 
and ® indicates the direction of the electric field vector E after transmission/reflection in 
accordance with our "rule of signs." 

The equations for refiection coefficients of s- and p-polarized waves from any complex 
multilayer (see Figure E]) also have a form similar to equation (j4]): 

'7into '7 

^ ^(1) -^(0) ,„. 

-^ 7into _|_ 7, ' V V 

where 2^(i)° is an apparent input impedance for a multilayer, i.e., it is the impedance that 
is seen by an incoming wave as it approaches to the interface. 

If the multilayer is made up of A^ plane-parallel, homogeneous, isotropic dielectric layers 
(with refractive indices Uj and geometrical thicknesses dj, where j = 1, 2, ... , A^) between 
semi- infinite incident (o) and external (e) media (see FigQ), the apparent input impedance 
Zy^\° of a semi-infinite external medium (e) and layers from A^ to j may be calculated by the 
following recursion relatioi^^^: 

„,„,„_„ Z;°S)-'%)tan(^) 

where aj = kz{j) dj = {2n/ X)nj cos{9j) dj] j = N,N — 1, ... ,2,1 and Zy^"^^^ = Z(^n+i) = ■Z'(e), 
riAT+i = He while dN+i = lie = by definition. 
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FIG. 3: Reflection and transmission for a multilayer with N layers. 



Fresnel's formulas for multilayer transmission coefficients are as follows: 
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where T(i+^\ are transmission coefficients at an interface between the j layer and the j ' + 1 
layer. 



B. The input impedance of a semi-infinite 1-D PC 

Now, we are ready to obtain the input impedance of the 1-D PC. Let us assume that we 
have a semi-infinite multilayer that consists of alternative layers with impedances Z(2) and 
Z(i) and its input impedance Z^(pc\ is unknown (see the top of Figure (|4])). To find it, we 
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FIG. 4: Determination of the impedance of the semi- infinite multilayer. Note that 
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make a next trick: first, we add an additional layer with impedance Z(i) to the multilayer 
and find the input impedance of this system using the recursion relation (j8]) : 
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Then, we add an additional layer with impedance Z(2) to this system and obtain the same 
semi-infinite multilayer again (see the bottom of Figure (jl])) with the input impedance: 



_ 7 ^(PCfc%))-^^(2)tan(a2) 
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It is obvious that Z^^p^^^ ^^^ ^.= ^'/p^x, and, by solving equations (IT2|) and (ITT]) for ZY^^y 

we obtain: 

^into ^ ((^?2) - ^fi)) tan(a2) tan(ai) ± y/j) 

(^^) 2 Z(2) tan(ai) + Z(i) tan(a2) ' 

where s = -4Z(i)Z(2) (^^(2) tan(ai) + Z(i) tan(a;2)j f^{i) tan(Q;i) + Z(2) tan(a2)) + 

(Z?2) ~ ^fi) ) tan(ai) tan(a2) • 

C. Dispersion relation for surface waves 

Now, we find the dispersion relation for surface waves in 1-D structures ended by an 
arbitrary layer with impedance Z(3). (It may be a layer with any na, for example a metal 
layer). The input impedance of such a structure will be: 

^into _7 Zjg°)-.Z(3)tan(a3) 

A general condition for the existence of a surface wave between two media with impedances 
Z\eit and bright is: 

■Z^left + bright = . (15) 

In our case, this condition takes the form: 

^(PC&Z(3)) + ^(e) = • (16) 

By solving equations (TTB1) and (TTll) . we obtain the dispersion relation for the optical surface 
waves in the 1-D PC: 



as 
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where M is a whole number. 

This is a general dispersion relation, which is valid for both polarizations. If one would 
like to obtain the dispersion of the s-polarized optical surface wave, one should use the Zs 
impedances in (TT7|1 and ( 1T^ . Accordingly, to obtain the dispersion of the p-polarized optical 
surface wave, one must use Zp impedances in (TT7|) and ( IT3l) . 

It should be noted that the solution obtained for s-polarization is equivalent (at M = 0) to 
the solution derived from a dispersion relation for an s-polarized optical surface wave, which 
was deduced by Yeh, Yariv and Hong using the unit cell matrix method (see equation (58) 
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in Ref.-). Meanwhile, the solution for p-polarization is equivalent to the solution derived 
from the dispersion relation for a p-polarized optical surface wave, which we deduced in 
a similar manner (using the unit cell matrix method) in Ref.- (see equation (2) in Ref.-). 
The dispersion relation flT7|) and the relations presented in Refs.-i^ were derived by different 
approaches (but, of course, both started from the same Maxwell equations and boundary 
conditions) and presented in different terms, but these relations give the same results and 
may be transformed to each other. The advantages of the presented dispersion relation with 
the impedance terms are its compact structure and its unified form for both polarizations. 
In addition, a visible physical interpretation the impedance terms in the current dispersion 
relation allows it to be easily extended for use with more complicated structures. For 
example, the addition of an adsorption layer between the layer (3) and the external medium (e) 
may be taken into account simply by changing of the impedance Z(e) by the impedance 
Z^^ calculated via recursion relation (|8]) , where impedance of the external medium Z(e) is 
convoluted with an impedance of the adsorption layer Z(a) ■ 

D. Band gap maximum extinction per length 

As a rule, in practical applications, we have values of two RIs of alternative media in the 
1-D PC, and the purpose is to find the thicknesses of each alternative layer, which provides 
the maximum extinction per length at given RIs, wavelength and angle. Below, we derive 
this condition for the maximum extinction, which makes it possible to minimize overall 
thickness of 1-D PC. Also, we show that the commonly-held opinion that it is a "quarter- 
wave-length" thickness of the layers that provide the maximum extinction per length is 
incorrect. 

To find the condition for the maximum extinction we derive the transmission coefficient 
for one period of the PC (two layers), that is, the transmission through the three interfaces: 

T/pcu+3)\ = T/pcu+i)\T(Pcu+2)\T/pcu+3)\ , (18) 

I PCU) ) \ PCU) ) \PCU + 1)) \PC(j+2)) 

where expressions for T/3+i\ are given by equations ([9]) or fITOl) . As a result we obtain for 
both polarizations: 



{Z(2) + ^|??S)) (^(1) + ^(pcfcZ(i))y .(<,,+„,) 



T.Pco.3)^ = f-^ )-'{ ' ' ^--^^^"^ e^(-+-) , (19) 



where Z^^pc^z i ) i^ given by (fTTj) and Z^(pq^ is given by (fT3|) . 

The desired values of the thicknesses d\ = (iimax and d2 = (i2max are the thicknesses at 
which the next expression: 



f{di,d2 



In 



\ PC(j) ) 



d\ + d2 



(20) 



reaches its maximum. Absolute values are inserted in equation (1201) in order to obtain the 
same result regardless of the ± sign in equation flT3l) . It may be noted that the "quarter- 
wave-length" thicknesses of the layers, dj{6j) = X/{4njCos{6j)), are values that maximize 
another expression, namely, expression fl2Ul) without the denominator, i.e., {di+d2)y<f{di, ^2) 
and these "quarter-wave-length" thicknesses are not optimal, especially at large incident 
(grazing) angles (i.e. at cos(6'j) — )■ 0). 

Now, we have all the equations required to calculate a 1-D PC structure for any particular 
experimental conditions. Presented below is an example in which we obtain the 1-D PC 
structure ended by a Pd nanolayer, which were successfully used for hydrogen detection in 
ReisM,^. 

III. CALCULATION OF THE 1-D PC STRUCTURE FOR PARTICULAR EX- 
PERIMENTAL CONDITIONS 

A. Angle-dependent variables 

It is very convenient to use a numerical aperture p = raosin(6'o) = ky{X/27T) as an angle 
variable instead of angles 6j in each j-th layer and hereafter we will do so. It is a unified 
angle variable for all layers, since according to Snell's law p = ?7,osin(6'o) = njsm{6j), for 
any j. The angle-dependent variables have the next forms as functions of p: 



Zs(4) = 77-T = , (for TE wave) (21) 



Z,(,)= ^^^ ^Vlzi^A!^ (for TM wave) (22) 



rij rij 



277 _ . 2tt 



kz{j) = —rij cos{6j) = —UjJl — {p/rijY (for both polarizations) . (23) 

A A 
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B. Conditions for long-range propagation of the plasmon-polariton waves in metal 
films 

A value of p at which we would like to excite a PC SW (and, therefore, a desired SW 
wavevector fcsw = ^irp/X) depends on the particular problem that we are trying to solve. As 
one example, to excite LRSPPs in thin metal films, the effective RI of the LRSPPs p should 
be close to the RI of the external medium Ue- To prove this statement, in this subsection 
we will find the angular value p at which the electric field of the incident p-polarized wave 
has a minimum inside a thin layer with a large extinction for optical waves. This minimum 
coincides with the zero of the main, tangential component Ey of the electric field in the 
film with the large extinction. This large extinction, i.e., large imaginary part of the RI of 
a thin film material {Im{nu) » 1), may arise both from a large negative permittivity of 
the material {Re^Su) « —1, even at small Im{eu)i as is the case for silver or gold) and 
from a large losses in the material {Im{£y[) » 1, even at small or positive Rei^EM), as is 
the case for palladium). In both cases the optical wave does not penetrate deeply in the 
bulk material with permittivity £3 = £m- It should be also noted that, in a very thin film, 
additional losses always appear due to collision-induced scattering of conducting electrons 
at the walls of the film. Therefore, the imaginary part of the permittivity in such a film is 
increased in comparison with the bulk material (see Appendix B for details). 

So, let us assume that we have a thin film with RI n^ = n-^ and that a p-polarized light 
wave is incident on it with angular parameter p. The instantaneous value of the tangential 
component of the electric field at a coordinate point z inside the film is the sum of the 
progressive and the recessive (reflected) waves: 

Ey{z) = ^(^)e-*'=-(3)('^«-^) + E(_)e*'=-(3)('^3-^) 

where the reflection coefficient from the interface between the film (3) and the external 
medium (e) is given by ([7]): 

To find the coordinate Zq at which the tangential component of the electric field in the 
film is zero, we solve the equation: 

Eyizo) = (26) 
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with respect to zq. Then we express the coordinate zq as zq = ds{l — a), where a is the 
coordinate of the zero minimum zq in terms of film thickness d^. When zq = (=> a = 1), 
the zero of the Ey takes place at an internal border of the film. When zq = d^/2 (=> a = 
1/2), the field zero occurs in the center of the film. When zq = d^ (=> a = 0), the field 
zero is located at an external border of the film. The solution of equations (126^12^ is: 



a = [1 - (zo/d^)] 



-i/fl(. 



2A;2(3)4 



(27) 



Assuming that |n3| >> p > rif, we will use the impedance of the final (metal) film 
in the so-called Leontovich approximation (see, for example, Ref.— for details). In this 
approximation equations fl22|) and fl23l) have the form: 
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(30) 



(I) 
p we get (in the approximation ^3 << A): 



and in equation f l27|) and then solving it for 
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Ue + 2nt 



an- 
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(31) 



which we presented earlier (see eq. (1) in Ref.-). 

So, we have derived equation (^11) . which shows that the modulus of the electric field 
strength \Ey\ has a minimum equal to zero when the p-polarized electromagnetic wave 
incident on the thin (metal) film in the angular range from pi to po- In this case, the point 
at which \Ey\ = is changed from the internal (a = 1, ^o = 0) to the external (a = 0, ^o = 1) 
borders of the film, and it is located in the center of the film at p = pi/2 (« = 1/2, zq = 1/2). 

Now the next question arises: are there some electromagnetic surface modes with a 
wavevector in the range: k = [-pi . . . -po]? If the answer is yes, one may expect that these 
modes will be long-range propagated surface modes, since they may be excited by the p- 
polarized waves incident on the film in the angular range p = [pi . . . po] , and these modes 
have the zero minimum of \Ey\ inside the lossy (metal) film. 
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One example of such a mode is widely recognized - it is the LRSPP in a thin film 
embedded between two identical dielectrics. It is known that the dispersion curve of SPPs 
splits as a result of a coupling between SPPs from both film interfaces, and the LRSPP 
wavevector shifts (at a given frequency) to a light curve (i.e., p — ?■ rie). The value of the 
LRSPP wavevector is fc ~ -pi/2 (see Ref.- and/or Appendix A), where pi/2 is given by (13 ip 
at a = 1/2, and the zero minimum of \Ey\ is always located in the center of the film. 

PC SWs are another example of electromagnetic surface modes that can propagate along 
thin metal films and have a wavevector in the range oi k = [-pi . . . -po\- In the next 
subsection we show how to design the PC structure for excitation PC SW at desired p ~ 
[pi . . . po] at given A, n\ and 77-2. 

C. Derivation of di, ^2 and da for PC structure terminated by Pd nanofilm 

So, let us to design PC structure with a maximum band gap extinction at desired p ~ 
Pi/2 — 1.0012 (obtained from fl3T|) at n^ = 1.0003, d^ ~ 10 nm and A ^ 739 nm). Let us 
have two layer materials, Ta20^ and Si02, with RIs of ^2 = 2.076 and rii = 1.455 at the 
wavelength of our tunable diode laser (at A ~ 739 nm). The function f{di, ^2) given by ( 120|1 
is presented in Fig. [5l One can see that f{di, ^2) reaches its maximum at dimax = 155.0 nm 
and (i2max = 112.8 nm. At these thicknesses, the maximum band gap extinction per length 
occurs, and, therefore, the PC structure will have a minimal overall thickness. 

Now, we must find the thickness d^ of the final palladium layer that will satisfy the 
dispersion relation ( TT7I) . With M = we have: 



2TTn^Jl - [p/n^^Y \ ^(3) + ^(PC)^(e) 



Substituting n^, = 1.9 + i4.8 and the other values pointed in the present subsection into f l32|) 
(and into (!22|) . (!23|) and (TT3|) correspondingly) we obtain ^3 = 1.2 nm. This is a rather small 
value and it is better for the thickness of the palladium layer to be in the range of 8 ... 10 nm 
to be sure that a continuous film we will obtain during deposition. If we want to retain the 
maximum band gap extinction at the desired pi/2 and, therefore, do not want to change 
the optimal values (iimax and (i2max, we may decrease the thickness of just the Ta20^ layer, 
which is contiguous with the Pd layer. The decrease in the thickness of the last Ta20^ layer 
from d2 = 112.8 nm to ^2 = 103.4 nm permits us to increase the thickness of the Pd layer 
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FIG. 5: 1-D PC band gap extinction as a function of layers thicknesses. 



to the value d^ = 8 nm and satisfy the SW excitation condition at the desired angle pi/2 
and the desired wavelength A. 

Thus we derive the next PC structure: substrate/{H L)^^ H' M /air, in which if is a Ta20^ 
layer with a thickness 62 = 112.8 nm, L is a Si02 layer with di = 155.0 nm, H' is a Ta205 
layer with cig = 103.4 nm, and M is the palladium layer with ^3 = dM = 8 nm. In Fig. [6] 
the calculated dispersion of this 1-D PC structure in air is presented as the logarithm of the 
optical field enhancement (i.e., as logT/eO in the external medium near the structure. The 
optical field enhancement Tfe\ is calculated by using ( ITOl) and (jH]) for 29 dielectric layers and 
one metal nanolayer on BK-7 {riQ = 1.513) substrate. The optimal number of Ta20^/ Si02 
pairs depends on a total extinction in the layers. In the present case, the Pd nanolayer is 
the layer that gives the main contribution in extinction and the 14 pairs provide the optimal 
coupling (i.e., Rpc = 0) of the incoming EM radiation into SWs. 

The dispersion is presented using the coordinate A(p). The angular parameter p, at 
which the excitation of the surface mode occurs, is equal to the effective RI of the mode, 
risw. Therefore, the red curve that is inside the blue band gap in Fig. |6]presents the dispersion 
of the LRSPP mode (i.e., the dependence of its effective RI on the wavelength). One can 
see from Fig. |6] that in accordance with fl3T|) the LRSPP mode exists near p r^ Uf. only 
(p :ii 1 . . . 1.0026), where the minimum of EM field inside the metal nanolayer occurs, while. 
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FIG. 6: The calculated dispersion of the 1-D PC structure with the terminal Pd nanolayer in air. 
The LRSPP mode is seen as the red curve (with an enhancement of more than 100) inside the 
band gap (blue areas with an enhancement of less than 1). The photonic band gap vanishes near 
Brewster's angle {pBr — 1-2 in this system), where no reflection of the TM wave takes place from 
the Si02/Ta205 interface. 

at p > 1.0026 the EM damping is increased significantly in this 8-nni thick metal film. 

In the presented case the optical surface mode dispersion curve approaches the line of 
the total internal reflection (TIR) (ptir = "^e = 1-0003) at wavelengths in the range of 
A ~ 734 . . . 746 nm. Therefore, at these wavelengths, we can excite the SPP at p — > rig 
and expect it to be LRSPP. At A = 733.7 nm (p = 1.0026 — pi), the zero minimum of 
\Ey\ occurs at the internal interface of the 8 nm palladium film, while at A = 740.2 nm 
(p = 1.00088 ~ P1/2) the zero minimum of \Ey\ occurs in the center of this 8-nm thick 
nanofilm. Moreover, at A = 745.6 nm (p = 1.000301, i.e., p — > po ~ Ue) we can excite 
LRSPP with \Ey\ = near the external interface of the film, and this wave will have the 
largest propagation length. Theoretical propagation lengths, i.e., propagation lengths that 
take into account the internal damping in the Pd nanofilm but do not take into account 
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the surface scattering of LRSPPs to photons, for these wavelengths are: L-j^^-j = 0.14 mm, 
-^740.2 = 0.32 mm and L745.6 = 5.4 mm. Therefore, contrary to intuitive expectations, 
modes with p = pq -^ n^ {\Ey\ = at the external interface of the film, i.e., modes with 
an asymmetric field distribution in the nanofilm) are more long-range propagated than the 
mode with p = pi/2 i\Ey\ = in the center of the film). More details about this issue are 
provided in Appendix A. 

IV. CONCLUSIONS 

The present article provides the theoretical background and the algorithm for the design 
of 1-D PC structures that support the propagation of optical surface waves. We have used 
the impedance approach, which permits calculation for s- and p-polarizations by the same 
equations. The main results of this work are the input impedance of semi- infinite 1-D PC (iT3ll 
and the dispersion relation of PC SWs ( TT71) . Equations ( TT9l) and ( l20l) are needed for the 
design of a multilayer structure with maximal band gap extinction at the minimal overall 
thickness of the structure. Equation (13T!) is important for understanding of the physical 
reasons for the propagation of the LRSPP even in thin films with large extinction, while the 
equations in Appendix A provide additional insight concerning why slightly asymmetrical 
structures provide more long-range propagation of the surface waves. 
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Appendix A: Dispersion relation for LRSPPs in the symmetrical configuration and 
increase of propagation length associated ^vith slight asymmetry in the system. 

The propagation of LRSPP in a symmetrical structure (metal film between identical 
dielectrics) has been considered in several publications-'-. Here we will consider this point 
based on the impedance approach and then will add slight asymmetry to the structure. 

The dispersion relation, A(p), for surface waves in a metal film between two semi- infinite 
dielectric media can be obtained easily, for example, from equation f lT7|) or equation f l32|) . in 
which the 1-D PC is replaced by a uniform dielectric medium with a RI tiq {Zy^^^ — )■ 2'(o)): 



For a pure symmetrical configuration [uq = rif, and Z(o) = 2'(e)], the general dispersion 
relation (lAip has the form shown below: 

-=t- f ^^Sffe I . (A2) 



^ 2714^3^1 - (p/n3)2 \ ^fs) + ^fe) 

By taking into account the Leontovich approximation, ( 128|) and ( l29l) . the equation can be 
simplified further: 

_ r^ arctan -j e \i^ e j f A3) 

A 2nd3n3 y n^^ + ris^Ue'^ - n^'^p^ j 

For a very thin film (at ^3 — ?■ 0) one can obtain two solutions for p from (IA3P : 



nl 



Plrspp - Ue . 



ds 



(A4) 



Tid^nl \n1 T{d-iRe{n\) 

^A TiRe{nl)d3 2\ 



PsRSPP - ^^ - ^r;^\:r - ^=F^ ■ (A5) 



The former equation presents an approximation for the long-range branch of the disper- 
sion relation of a symmetrical system flA2p . while the latter is an approximation for the 
short-range branch (its damping increases at c?3 — )■ 0). In Fig. [71 graphical illustrations 
of dispersion relation ( JA2I) [or ( JA3I) . since no difference can be seen between them at this 
scale] and its approximations (1A4I) and ( 1A5I) are presented. Calculations are performed for 
an 8-nm thick Pd film that is suspended freely in the air (rig = 1.0003), and the Pd RIs 
at different wavelengths were taken from Ref.— . In addition to good agreement between 
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the dispersion equation for the symmetrical system and the approximated equations, (lA4p 
and flA5l) . Fig. [7] shows that the long-range branch of the curve has no noticeable wavelength 
dispersion in the visible range, while the system with PC from one side the LRSPP curve 
(see Fig. [H]) has prominent dispersion. 
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FIG. 7: The dispersion of a Pd nanofilm that is freely suspended in the air. The solid black 
curve is calculated from the general dispersion relation for the symmetric system (lA2h . while the 
dotted red curve and the dashed-dotted blue curve are calculated by approximated solutions ()A4p 
and (jASp . respectively. 



It is apparent that flA4p consides exactly with (!3T|) at a = 1/2, and, therefore, in a pure 
symmetrical configuration, the zero of the tangential component and the total minimum 
of the electric field are always located in the center of the metal nanofilm, which is not 
surprising due to the symmetry of the system. A more surprising fact is that this symmetrical 
system (with l^'yl = in the center of the film) is not the optimal solution if one is looking 
for the smallest damping and the longest propagation length of the LRSPP. To the best 
of our knowledge Wendler and Haupt were the first to point out on this peculiarity. In 
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the theoretical article^^ they showed, through numerical calculations, that, in a slightly 
asymmetrical structure, the propagation length may be three orders of magnitude greater 
than the propagation length in the symmetrical structure. 

In experimental practiceS,"i^ we have been able to increase the propagation length of 
SPPs by a factor in the range of 100-200 times. Further increase is limited by damping 
due to surface scattering, which occurs due to the close proximity of the effective RI of the 
LRSPP to the RI of the external medium, or, in other words, "the plasmon curve becomes 
too close to the light line". Therefore, any small disturbance can transform plasmons to 
photons when this zero minimum is approaching the border of the film (and, therefore, 
p — )■ rie). But, nevertheless, the internal damping of LRSPPs can be decreased when slight 
asymmetry exists in the structure. 

Below, we obtain analytical expressions for the optimal Anmax difference between the RIs 
of the dielectrics and show that the propagation length of the LRSPP increases when the 
minimum of the electric field shifts to the interface of the film (but does not leave the borders 
of the film). To do this, we return to the general dispersion relation for a non-symmetrical 
system ( JAll) and assume that the RIs of the dielectrics on both sides of the metal film differ 
only by a small value An: 

riQ = ne + An . (A6) 

We are solving for the LRSPP wave only, so p will differ from uq by a small value Ap: 

p = no + Ap = ne -h An -F Ap , (A7) 

where Ap ~ 0{dl/\^) - see f lAi]) . 

In the limit, with small Ap and small An, the dispersion relation (lAll) can be simplified 
to the form: 



A-^(Ap,An) 



' ^3 ^Ap+ -An) -{^nl + ^nl] ApAn 



(A8) 



Solving this equation for Ap and taking its imaginary part we obtain the next extinction of 
LRSPPs in a slightly asymmetrical system: 

Ini [pLRSPp(An)] = Im [no + Ap] ^ 

^J^}!^^^ (iQTi'^dlnl + Svr^rf^An (SRe{nl) + 9n^) - 8AnA^) . (A9) 

32nj^A^ ^ v / / 

This over-approximated imaginary part of the effective RI of the LRSPP has two solutions 
with (approximately) zero extinction, /m [pLRSpp(An)] = 0, and, therefore, with maximal 
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propagation length, namely, the first solution is An = and the second solution is: 

^'''""^ "^ ^^8A2-37r2cii(8L'(ni) + 9n2) ' ^^^^^ 



which in the limit rf'^ ?^ is: 



An^axo ~ 2n^ 



(ia 



'3 



(All) 



One can see that flAlip coincides with the addition to Ue from fl3T|) at a = 1 (i.e., when the 
minimum of the field is at the nanofilm border). Therefore, expression (JAIOP also may be 
considered as a cutoff condition for the existence of LRSPPs (at An > Arimax no bounded 
modes exist). From (lAip . it can be determined that, at rie = tiq — Arimax, the effective RI 
of LRSPP is infinitesimally close to the TIR angle, i.e., p — )■ uq. 

It is worth noting here that the analytical expression for Arimax from flAlOp corresponds 
well with numerical data presented by Wendler and Haupt in Ref.— . Even the over- 



approximated expression for Arimaxo from flAlip . which is not dependent on the nanofilm 



RI n3, is coincident with these numerical data at small film thicknesses (i.e., (is < 15 nm). 

At excitation of ultra-long-range SPP in thin films by introducing a small asymmetry 
in the system and, therefore, by abutting the plasmon curve on the light line, it is very 
important to maintain an appropriative balance between decreasing the internal damping 
of LRSPPs and increasing the scattering of the LRSPPs. The possibility for tuning system 
parameters during the experiment to find the optimum is extremely convenient in this case. 

The principal difference between D/M/D' and D/M/PC systems is that, in the system 
with a 1-D PC on one side, it is experimentally easy to introduce such a small asymmetry 
(in effective RI) by tuning the wavelength of the laser, while it is really difficult to (finely) 
tune the dielectric constants on the one interface by any means. Wavelength tuning would 
be not useful in the D/M/D' case, because dispersion of dielectrics is small in non-absorbing 
wavelength regions, while wavelength tuning in the PC structure is very effective due to the 
high wavelength dispersion inside the band gap regions. 
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Appendix B: Additional damping in thin film. 

Damping of EM waves in thin metal films increases due to collision-induced scattering of 
conducting electrons at the walls of the film. This occurs when the mean free path of the 
electrons becomes comparable to the characteristic dimensions of the system. Therefore, 
the imaginary part of the permittivity in such a film is increased in comparison with the 
bulk material. Unfortunately in many articles (especially theoretical articles), authors have 
usually disregarded this fact and used overly optimistic values for the imaginary part of the 
metal permittivity. 

Here we present a summary of formulae for this case, which help to estimate the minimal 
addition to Im^eu) = e'^^^^^k- ^^ optical frequencies, the damping in nanostructures with the 
characteristic dimension L is 

7 = 7buik + -^ , (Bl) 

where 7buik is the damping constant for the bulk sample and vp is the electron velocity on 
the Fermy surface. Below, several specific cases are considered. 

1. Nanofilm as a set of nanoparticles: 

For sphere^^ 

^ = ^7 = ^ (,7buik + ^)= .,.,, + -^-,^ , (B2) 

where Up is plasma frequency. Therefore, the characteristic dimension in this case is L = 
Lji = 4_R/3 {R is the radius of the sphere). 

2. Continuous film: 

For continuous film with thickness dr^: 

SujpVp 1 + cosh'^ (upd/c) 

7 = 7buik + -5 r-TT — TT^; TT — ,/ x , JT ■ 1^3) 

8 c sinh[ujpa/c) cosh[ujpa/c) + Upd/c 

For very thin film (nanofilm), in the limit (i — )• 



3 vp ,^ . . 

7 = 7buik+gY' ^^^) 
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(therefore L = L^ = 8d/3 in this hmit) and 



_„ Up (_ ^ 'iVF\ _„ , ?>UjIvf 



^-^l^- + WJ-^-+8^f' (««' 



while in the general case 

e" = e" + M!^ 1 + cosh'^ jupd/c) 

" 8u^ c sinh{upd/c) cosh.{upd/c) + Upd/c 

We have considered only collision-induced damping, and we did not take into account 
other types of damping in small volumes, which may lead to further increases of the imagi- 
nary part of the nanosized materials at optical frequencies. 

The collision-induced damping increase in the nanofilm at zero frequency may be es- 
timated by Fuchs' formula-^ for a specific electrical resistance of continuous films. For 
example, in our work with gold nanofilm^, the electrical resistance of the nanofilm was 850f2 
for a strip that was with 0.25 mm wide and 6.5 mm long. Therefore, for our nanofilm 
thickness of 5 nm, the specific resistance in the nanofilm was pfum = IG.S/ificm, which is 
7.4 times greater than the specific resistance of bulk gold (2.2/if2cm). While according to 
Fuchs formula at d << l^ 

Pj^ = ^xJh\, (B7) 

Pfilm 4/o V " / 

Substituting the mean free path of electrons in bulk gold, which at room temperature is 
/o — 50 nm, we obtain that the specific resistance of the 5 nm gold film at room temperature 
should be about six times greater than that of bulk gold. This value is in good agreement 
with the measured value, taking into account imperfections of the sputtered gold nanofilm. 
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